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•  Efficient  Composite  Designs  with  Small  Number  of  Runs 
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0  Summary 

A  new  method  of  constructing  composite  designs  using  the  robustness  property  against 
deletion  of  runs  (points)  is  given.  Composite  designs  are  presented  with  small  number  of 
runs.  These  designs  are  more  efficient  in  terms  of  both  prediction  of  response  and  estimation 
of  parameters,  than  their  competitors  that  are  available  in  the  literature. 

Short  Running  Title:  Composite  Designs 

Key  Words:  Deletion  of  Points,  Factorial  Points,  Orthogonal,  Response  Surface,  Robustness, 
Unbiasedness. 

AMS  Subject  Classifications:  Primary  and  Secondary  62K15,  62J05 

1  Introduction 

In  the  study  of  dependence  of  a  response  variable  y  on  k  explanatory  variables,  coded 
as  Xi,...,Xk,  the  unknown  response  surface  is  approximated  by  a  first  or  a  second  order 
polynomial  in  a  small  region  with  the  center  being  the  point  of  maximum  interest.  If  the 
first  order  model  gives  a  significant  lack  of  fit  indicating  the  presence  of  a  surface  curvature,  a 
second  order  response  surface  model  is  then  fitted.  The  N  runs  (or  points)  in  the  design  are 
(xui, . . . ,  xuk)  and  the  observations  are  y(xu i, . . . ,  xuk )  =  yu,  u  =  1, . . . ,  N.  The  expectation 
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of  yu  under  the  second  order  model  is 

k  k  k  k 

E(V u)  =  fio  +  J2  /%*«<  +  £  P»xli  +  £  PijXuiXuj,  (1) 

t=l  i=l  1=1  j=l 

«<  i 

where  the  intercept  fi0,  the  linear  coefficients  the  pure  quadratic  coefficients  fin  and  the 
interaction  coefficients  fiij  are  unknown  constants.  The  yu’s  are  assumed  to  be  uncorrelated 
with  the  variance  a2,  an  unknown  constant.  The  number  of  fi's  in  (1)  is  1+2 k  + 

For  a  second  order  design  with  N  points  (xui? •  •  •  ,  xufc)iu  —  1, ...  ,N,  all  fi' s  are  unbiasedly 
estimable. 

A  special  second  order  design,  called  composite  design  (CD),  consists  of  F  factorial  points 
(FP’s)  which  are  a  fraction  of  2k  points  (±1, . . . ,  ±1),  2k  axial  points  (AP’s)  (±a, . . . ,  0), . . . , 
(0,...,±a),  a  is  a  given  constant  and  n0(>  0)  center  points  (CP’s)  (0,  ...,0).  The  total 
number  of  points  is  N  =  F  +  2k  +  n0.  Box  and  Wilson  (1951)  introduced  such  designs. 
Box  and  Hunter  (1957)  suggested  FP’s  as  the  complete  set  of  2k  points  or  an  orthogonal 
resolution  V  plan  (i.e,  the  plan  that  permits  the  unbiased  estimation  of  fi0  +  fin-,  fii's 
and  fiifi s  under  (1)  and,  moreover,  the  estimators  are  un correlated).  The  CD’s  with  such 
FP’s  give  the  variance  of  the  predicted  response  dependent  on  the  point  only  through  its 
distance  from  the  origin.  This  variance  structure  is  achieved  at  the  cost  of  a  large  number 
of  points,  particularly  FP’s.  Efforts  are  then  being  made  for  reducing  the  number  of  FP’s. 
Hartley  (1959)  pointed  out  that  FP’s  need  not  be  of  resolution  V  but  could  be  as  low  as  of 
resolution  III  plan  (i.e.,  the  plan  that  permits  the  unbiased  estimation  of  fiQ  +  fin  and 
fii's  assuming  /?,/ s  are  known)  with  an  additional  condition  that  the  unbiased  estimation 
of  fi^' s  is  possible  assuming  the  other  fi' s  are  known.  Draper  and  Lin  (1990)  named  such 
FP’s  as  resolution  IIP.  Hartley  (1959)  presented  resolution  IIP  FP’s  as  regular  fractions 
and  Westlake  (1965)  presented  FP’s  as  irregular  fractions  of  2k  factorials.  Draper  (1985), 
Draper  and  Lin  (1990)  gave  resolution  IIP  FP’s  using  the  projection  properties  of  Plackett 
and  Burman  orthogonal  resolution  III  fractions  of  2k  factorials. 
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In  this  paper,  a  new  method  of  constructing  CD’s  is  given  by  introducing  first  a  sub¬ 
model  of  (1),  presenting  orthogonal  FP’s  under  the  submodel  and  finally,  reducing  the  num¬ 
ber  of  FP’s  using  the  idea  of  robustness  of  designs  against  deletion  of  points  [see  Ghosh 
(1979)].  New  CD’s  constructed  by  this  method,  can  be  made  minimal  in  the  number  of 
FP’s.  They  are  then  compared  with  the  available  CD’s  in  the  literature.  Comparisons 

i 

are  made  with  respect  to  the  Trace,  Determinant  and  Maximum  Characteristic  Root  of  the 
variance-covariance  matrix  of  the  least  squares  estimators  of  /?’ s  and  also  with  respect  to 
the  average  of  the  variances  of  predicted  responses  in  a  spherical  region  about  the  center. 
Our  designs  perform  significantly  better  over  the  designs  that  are  available  in  the  literature. 
Results  presented  in  this  paper  are  striking,  useful  and  valuable. 

2  Factorial  Points 

Observations  at  AP’s  permit  the  unbiased  estimation  (UE)  of  /?,•  and  f30  +  a2 i  = 
1, . . . ,  k.  The  CP  observations  provide  the  UE  of  /?„.  For  a  second  order  CD,  observations 
at  FP’s  must  at  least  allow  the  UE  of  /?0  +  J2i=i  At  and  /?,j,  i  <  j,i,j  =  1  ,...,&,  given 
that  the  estimators  of  are  available  from  AP’s.  When  there  is  no  center  point 

observation  (i.e.,  na  =  0),  a  can  not  be  equal  to  k1^2.  In  view  of  this,  the  following  submodel 
of  (1)  is  introduced  for  the  choice  of  FP’s. 

k  k 

E{Vu)  =  &  +  ]££ PijXuiXuj,  U  =  1,...,F.  (2) 

i= 1  j= 1 
*<j 

In  matrix  notation, 

E(y)  =  X&  V(y)  =  o*T,  (3) 

where  y(F  x  1)  is  the  vector  of  observations  at  FP’s,  ^(p  x  l),p  =  1  +  Q),  is  the  vector  of 
/3’s  in  (2)  and  X(F  x  p)  is  the  design  matrix  based  on  FP’s.  It  is  important  to  note  that  a 
CD  with  a  >  0  for  n0  >  0  and  in  addition,  a  ^  &1/2  for  n0  =  0  is  of  second  order  if  and  only 
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if  the  UE  of  all  /?’ s  under  (2)  is  possible  (i.e.,  Rank  X  =  p ). 


Let  Fw  be  a  set  of  FP’s  with  the  corresponding  vector  of  observations  y  ,  w  =  1, 2,  and 
£  =  where  £(p,  xl),  i  =  1,2,3,  with  p, +pj +p,  =  p.  Denote 


where  Xwi(Fw  x  p,-),u>  =  1,2,  z  =  1,2,3.  Let  F\  and  F2  FP’s  be  such  that 

X12  =  Xiz,X22  =  —X23,  P2  =  P3- 


(5) 


Example  1.  For  k  =  8,  the  64  FP’s  satisfying  x4x2x3  =  x4x5x6  =  -1  form  OFP’s.  For 
F\  —  32  FP’s  satisfying  xi£2X3  =  x4x^Xq  =  —1,  x4x4x7x&  =  1  and  F2  —  32  FP’s  satisfying 
X4X2X3  =  *4X5*6  =  x4x4x7x6  =  -1,  g  =  (014,017,04?)  and  =  (07S,04 8,/318),  p2  =  p3  = 
3,  pi  =  23,  the  conditions  in  (5)  in  fact  hold.  The  following  result  is  very  useful  in  the 
selection  of  FP’s. 

Theorem  1.  If  for  F1  FP’s  Rank  [*n,*i2]  =pi+p2  and  for  F2  FP’s  Rank  X22  =  p2  then 
for  (Fx  +  F2)  FP’s, 


Rank 


- 1 

* 

M 

*12 

1 - 

CO 

[*21 

*22 

* 

00 

1 _ 

=  p. 


Proof.  There  exists  a  p2  x  p2  submatrix  X^2  of  X22  with  Rank  X£2  =  p2.  The  matrix 
obtained  by  taking  the  corresponding  p2  rows  of  [X2U  X22,  X23]  is  [X^,  X^X^].  It  now 
remains  to  prove  that 


Rank 


Xu  X42  X\3 


Y *  Y*  Y* 
L'a21  ^22  A23j 


=  p. 


Let  if 
matri 


A  set  of  FP’s  is  said  to  be  orthogonal  if  X'X  is  a  diagonal  matrix  and  is  denoted  by 
OFP’s.  FP’s  which  are  not  OFP’s  are  called  nonorthogonal  FP’s  (NOFP’s).  Notice  that  if 
one  FP  is  negative  of  another  FP,  then  the  corresponding  rows  of  X  are  identical. 

Thus 
from 
D2  = 
(4)  Rani 


Exar 

Rani 
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Let  if  possible  that  rank  be  (p— s),  where  s  >  0.  Thereexists  an  ((F1  +  F2  —  p+-s)x(i7\+ii2)) 
matrix  [D\,  D2]  with  rank  equals  to  F\  +  jF2  —  p  +  s  and  satisfying 


[DuDt] 


Thus  DiXu  +  D2X22  —  0  and,  moreover,  D^X^  +  Z>2A£3  =  DtXu  -  D2X22  -  0  (since, 
from  (5),  Xu  -  X13  and  X\2  =  -X^).  This  implies  that  D2X^  =  0  and  consequently 
D2  =  0.  Then  Rank  Di  =  (Fi  +  F2  -  p  +  s)  >  (Fa  -  pi  -  p2)  and  this  is  impossible  since 


Rank  [Xn,Xi2  ,  X13]  =  pi  +  P2-  This  completes  the  proof. 


Example  2.  For  k  =  8  and  Fi  =  32  FP’s  satisfying  xix2x3  =  x4x5x6  =  -1  and  xix4x7x8  =  1, 
Rank  [^11,^12]  =  Pi  +  P2  —  23  +  3  =  26.  From  (5),  X22  =  — X23  and  the  columns  of  X22, 
correspond  to  /?18,  /?48  and  /?78.  It  is  now  clear  that  the  condition  Rank  X22  —  P2  =  3  in 
Theorem  1  can  be  achieved  by  choosing  3  points  from  32  points  satisfying  x12r2®3  =  X4X5X6  = 
x\x4x7xs  =  —1  so  that  the  columns  for  x\,x4  and  x7  are  independent.  For  illustration,  one 
such  choice  for  F2  =  3  FP’s  is 

--1  1  1  1-11  1  1- 

11-1-1  11  11. 

II;  1  1  — 1  1-11-1  1. 


A  general  series  of  F\  and  F2  FP’s  satisfying  (5)  is  given  below. 


I.  fc  =  3i  —  1,  t  >  2  and  t  is  an  integer. 

The  Fu,  u  —  1,2,  FP’s  satisfying 

XjX2X3  =  ...  —  2^3t— 5X3t— 4X31— 3  t  ^  2, 

2:1X4X73:8  =  •  •  •  =  a:t+3[t/3]_5Xi+3[1/3]_2X3t-23;3i-l  =  (3  -  2 u),t  >  3, 


k  =  4,  x4  =  — 1,  F  =  8. 

/:  =  5,  Plan  I.  2122X324X5  =  —1, 

Plan  II.  x5  =  —1, 

F  =  16  for  both  plans. 

k  =  6,  21X2X3X4X5  =  x6  =  —  1,F  =  16.  (8) 

The  (Fi  +  F2)  FP’s  in  (6)  and  (7)  and  the  F  FP’s  in  (8)  form  OFP’s  under  (3).  It  follows 
from  (6)  that  the  F\  and  F2  FP’s  for  k  =  St  +  1,  are  in  fact  twice  of  their  counterparts  for 
k  =  3 1.  The  Fu ,  u  =  1,2,  FP’s  for  k  =  3t  +  1,  can  be  obtained  from  the  corresponding  FP’s 
for  k  —  3t  by  adding  1  and  —1,  respectively  in  the  last  position. 


str 
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3  Construction  of  Minimal  Plans 


The  number  of  FP’s  discussed  in  Section  2  is  in  fact  large.  The  idea  of  robustness  of 
pp’s  against  deletion  of  points  [Ghosh  (1979)]  is  now  used  in  reducing  the  number  of  FP’s. 
The  idea  turns  out  to  be  extremely  powerful. 

Dpfinition.  A  set  of  F  FP’s  is  said  to  be  robust  against  deletion  of  t  (a  positive  integer) 
points  if  the  parameters  in  /?  under  (3)  are  still  unbiasedly  estimable  with  the  remaining 
(F  - t)  FP’s. 

For  a  robust  set  of  F  FP’s,  the  rank  of  the  resulting  matrix  remains  p  when  t  rows  of  X 
corresponding  to  t  points  are  deleted.  The  following  result  of  Ghosh  (1979)  is  instrumental 
in  determining  the  robustness  of  FP’s. 

Theorem  2.  Let  Z((F  -  p)  x  F)  be  a  matrix  with  Rank  Z  =  (F  -  p)  and  ZX  =  0.  A 
necessary  and  sufficient  condition  for  a  set  of  F  FP’s  to  be  robust  against  deletion  of  t  points 
is  that  Rank  Z*  =  t ,  where  Z*{{F-p)  x  t)  is  the  submatrix  of  Z  corresponding  to  t  points. 

The  following  corollary  is  very  useful  in  determining  robustness  by  exploiting  the  design 
structure. 

Corollary.  Let  X'  =  [X^X^]  where  Xr(fi  x  p),  X2{f2  x  p),  fi  +  h  =  F  and  ^1^2  =  °- 
If  the  set  of  /,  FP’s  is  robust  against  deletion  of  U,i  =  1,2,  points,  then  the  set  of  F  FP’s  is 
robust  against  deletion  of  (t\  +  £2)  points  and  vice  versa. 

Proof.  Note  that  Rank  X  =  Rank  Xx  +  Rank  X2  and  the  rest  is  clear. 

Plans  are  now  constructed  for  4  <  k  <  10. 
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3.1  (k  =  4). 


Two  plans  with  8  PF’s  in  (6)  and  (8)  are  robust  against  deletion  of  any  one  point.  The 
resulting  plans  with  7  FP’s  are  denoted  by  Plan  4.1  and  Plan  4.2.  There  are  8  possible 
plans  for  each  of  Plan  4.1  and  Plan  4.2. 

i 

3.2  (k  =  5). 

Consider  F1  -  16  FP’s  given  in  (6).  Note  that  p2  =  p3  =  0.  The  rows  of  the  matrix 
Z(5  x  16)  in  Theorem  2  correspond  to  x4,  x5,  X1X4X5,  x2x4x5,  x3x4x5  and  the  columns 
correspond  to  16  points.  It  follows  from  the  structure  of  Z  that  the  5  points  with  exactly 
one  of  the  x,-’s  equal  to  1  or  —1,  can  be  deleted.  The  resulting  11  points  satisfy  Rank  X  =  11 
under  (3). 

Consider  F  =  16  FP’s  given  in  (8)  satisfying  x4x2x3x4xh  =  —1.  The  11  points  obtained 
by  deleting  the  5  points  with  exactly  one  of  x.-’s  equal  to  -1,  satisfy  Rank  X  =  11. 

Similarly  consider  F  =  16  FP’s  given  in  (8)  satisfying  **  =  -1.  The  11  points  obtained 
by  deleting  the  5  points  with  exactly  one  of  x.’s  equal  to  1  or  -1,  also  satisfy  Rank  X  =  11. 
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Another  set  of  11  FP’s  is  also  given  below. 

-1  -1  -1  -1  -r 
1  -1  -1  -1  1 
1-1  1  1-1 

1  1-1  1-1 

1-1-1  1  -1 

1-1  1-1  1  .  (9) 

1  1-1-1  1 

111-1-1 
1  1-1-1  -1 

1-1  1-1  -1 

-1  1  1-1  -1. 

The  above  4  plans  with  11  FP’s  give  identical  X  matrix  under  (3).  The  plan  given  in 
(9)  is  named  as  Plan  5.1. 

3.3  (k  =  6). 

The  plans  with  16  FP’s  given  in  (6)  and  (8),  are  called  Plan  6.1  and  Plan  6.2,  respectively. 
These  plans  have  the  minimum  number  of  points  under  (3). 
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Table  1 


Efficiency  Comparison  for  k  =  4. 


nQ 

a 

Plan 

MCR 

T 

D 

V{r) 
r  =  1 

matm 

■ 

V(T) 

r  —  \fk 

0 

i 

BB 

10.194 

1.668e-09 

23.194 

70.139 

7.194 

1.657e-09 

19.444 

55.139 

1 

i 

4.589 

9.843 

1.287e-09 

10.345 

24.131 

73.262 

gfUl 

1.448 

7.035 

1.277e-09 

9.409 

20.387 

58.284 

m 

4.1 

13.151 

19.417 

m 

4.2 

13.089 

18.417 

4.1 

4.414 

9.626 

1.052e-09 

10.535 

25.237 

76.818 

4.2 

1.436 

6.935 

1.038e-09 

9.582 

21.427 

61.576 

mm 

0.988 

8.602e-16 

8.041 

9.673 

■ 

0.783 

2.898 

9.081e-16 

7.974 

9.408 

18.151 

9.214 

8.151 


riQ  ct  Plan 


Table  2 

Efficiency  Comparison  for 

MCR I  T  I  D 


k  =  5. 
V(r) 
r  =  1 


12.569 

15.687 

13.812 

13.153 

86.340 

99.125 

63.743 

67.528 

12.752 

15.967 

14.121 

13.418 


V(r) 
r  =  \/2 

V{r) 

r  =  sfk 

27.215 

107.236 

39.664 

185.937 

32.263 

'135.896 

29.611 

119.819 

62.111 

40.007 

73.250 

64.625 

49.056 

45.410 

50.861 

41.694 

28.247 

112.036 

41.150 

193.559 

33.570 

142.321 

30.785 

125.443 

Table  3 

Efficiency  Comparison  for 


MCR  T  D 


1.032 

1.032 

49.889 

78.443 

1.032 

1.032 


.168 

.168 


1.149 

1.164 

1.032 

0.493 

0.584 

0.585 

0.582 

0.587 


T 

D 

9.587 

9.227 

4.440e-25 

4.397e-25 

51.966 

80.402 

1.021e-34 

1.060e-34 

9.571 

9.191 

3.963e-25 

3.837e-25 

3.181 

3.069 

9.466e-37 

6.196e-37 

3.226 

3.123 

2.352e-36 

1.574e-36 

9.558 

2.917 

3.576e-25 

1.591e-34 

2.597 

2.486 

4.733e-37 

3.103e-37 

2.659 

2.546 

1.191e-36 

7.938e-37 

k  =  6. 


V(r) 
r  =  1 


16.211 

16.092 

907.747 

414.564 


16 

16 


V(r) 

r  =  V 2 
35.444 
34.968 


24.102 

24.066 

16.83 

11.473 

13.984 

13.848 

13.905 

13.768 


22.077 

21.756 

37.593 

13.781 

14.686 

14.368 

14.758 

14.393 


V(r) 

r  =  yfk 
176.900 
178.216 


42.893 

98.282 


34.968 

33.523 


36.018 

34.458 

189.482 

43.875 

34.674 

33.179 

35.912 

34.326 


no  a 


Table  4 

Efficiency  Comparison  for 
MCR  I  T  I  D 


7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 

7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 

7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 

7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 


7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 


7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 


7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 


7.1 

7.2 
7.W 
7.DL.1 
7.DL.2 


5.512 

9.246 

12.946 

58.800 

24.485 

13.232 

8.486 

27.740 

15.293 

10.871 

5.461 

9.246 

12.814 

58.800 

24.433 

2.211 

2.351 

1.926 

16.170 

7.628 

2.208 

2.211 

1.879 

15.286 

7.363 


5.419 

9.246 

12.708 

58.800 

24.391 


2.211 

2.351 

1.926 

16.170 

7.628 


2.186 

2.211 

1.804 

15.286 

7.351 


16.298 

20.184 

29.464 

77.428 

40.931 

17.672 

13.296 

33.611 

26.865 

21.521 

16.230 

20.173 

29.307 

77.418 

40.864 


5.790 

6.224 

7.469 

21.684 

12.736 

5.540 

5.810 

7.081 

20.390 

12.121 


16.178 

20.164 

29.182 

77.409 

40.807 


5.217 

5.652 

6.897 

21.113 

12.165 


5.035 

5.343 

6.528 

19.935 

11.623 


4.907e-33 

4.692e-33 

3.381e-31 

2.741e-31 

2.022e-31 

1.470e-48 

1.395e-48 

1.041e-46 

8.256e-47 

6.534e-47 

4.265e-33 

4.269e-33 

3.025e-31 

2.505e-31 

1.842e-31 


l.llle-48 

1.787e-48 

4.443e-47 

1.127e-46 

6.905e-47 


1.139e-49 

1.640e-49 

4.193e-48 

1.084e-47 

6.581e-48 


3.969e-33 

3.907e-33 

2.751e-31 

2.292e-31 

1.660e-31 


5.466e-49 

8.937e-49 

2.221e-47 

5.644e-43 

3.456e-47 


5.941e-50 

8.731e-50 

2.139e-48 

5.806e-48 

3.467e-48 


k  =  7. 
V(r) 
r  =  1 
22.472 
23.886 
27.260 
44.701 
31.430 
337.852 
220.443 
691.102 
199.638 
264.983 
22.785 
24.296 
27.668 
45.689 
32.018 


V(r) 

r  =  y/2 
53.007 
58.608 
72.179 
141.865 
88.799 

273.909 
180.273 

558.909 
180.609 
223.540 
54.212 
60.047 
73.798 
145.615 

90.994 


V(r) 

T  =  y/k 

387.556 

448.182 

625.412 

1467.853 

820.326 

137.261 

120.216 

238.511 

375.441 

244.855 

397.247 

460.609 

639.998 

1508.606 

842.064 


103.159 

126.286 

386.177 

222.437 


23.141 

55.454 

24.735 

61.508 

28.125 

75.472 

46.706 

149.388 

32.640 

93.221 

18.679 

22.148 

18.934 

22.956 

19.371 

24.817 

24.866 

46.684 

21.401 

32.903 

18.072 

21.535 

17.878 

21.771 

19.076 

24.274 

23.331 

44.025 

20.438 

31.530 

407.083 

473.040 

654.981 

1549.362 

863.930 


96.340 

104.411 

128.164 

395.079 

226.913 

92.825 

98.628 

121.153 

372.918 

216.696 


Table  5 


Efficiency  Comparison  for  k==_8. 


Plan 

MCR 

T 

D 

V(r) 
r  =  1 

V(r) 

r  —  \/2 

V(r) 
r  =  \Z% 

8.1 

4.507 

20.696 

1 .51 le-44 

28.511 

68.598 

640.190 

8.2 

4.454 

21.007 

1.436e-44 

28.633 

68.993 

635.767 

8.3 

4.456 

21.043 

1.707e-44 

28.645 

69.040 

636.517 

8.4 

4.455 

21.038 

1.419e-44 

28.645 

69.040 

636.517 

8.DL.1 

9.047 

29.812 

4.568e-43 

31.934 

82.196 

846.954 

8.1 

4.505 

20.648 

1.381e-44 

28.869 

69.885 

653.340 

8.2 

4.454 

21.000 

1.336e-44 

29.038 

70.349 

649.878 

8.3 

4.456 

21.035 

1.568e-44 

29.050 

70.396 

650.642 

8.4 

4.455 

21.032 

1.326e-44 

29.050 

70.397 

650.644 

8.DL.1 

9.047 

29.802 

4.223e-43 

32.411 

83.841 

865.722 

8.1 

1.212 

6.047 

7.592e-64 

41.226 

40.897 

132.250 

8.2 

1.141 

6.286 

1.191e-63 

41.400 

41.400 

138.000 

8.3 

1.126 

6.260 

1.074e-63 

41.391 

41.364 

137.425 

8.4 

1.136 

6.291 

l.llle-63 

41.403 

41.412 

138.192 

8.DL.1 

1.956 

8.756 

4.616e-62 

42.338 

45.166 

198.415 

8.1 

4.504 

20.608 

1.272e-44 

29.259 

71.200 

666.615 

8.2 

4.454 

20.994 

1.242e-44 

29.467 

71.722 

663.992 

8.3 

4.456 

21.029 

1.462e-44 

29.479 

71.771 

664.769 

8.4 

4.455 

21.025 

1.227e-44 

29.479 

71.771 

664.772 

8.DL.1 

9.046 

29.793 

3.885e-43 

32.912 

85.505 

884.496 

8.1 

1.212 

5.485 

3.802e-64 

23.078 

26.707 

133.558 

8.2 

1.141 

5.724 

5.962e-64 

23.255 

27.221 

139.433 

8.3 

1.089 

5.698 

5.379e-64 

23.246 

27.184 

138.846 

8.4 

1.135 

5.728 

5.500e-64 

23.258 

27.233 

139.629 

8.DL.1 

1.956 

8.194 

2.310e-62 

24.214 

31.069 

201.162 

Table  6 


Efficiency  Comparison  for  Jc  =  9. 


no 

a 

Plan 

MCR 

T 

D 

h-1  N 

V(r) 
r  —  72 

F(r) 
r  =  7k 

0 

i 

9.1 

1.873 

17.967 

2.486e-58 

32.383 

74.351 

746.289 

9.2 

2.755 

18.965 

2.140e-58 

32.767 

75.889 

777.443 

9.3 

3.946 

20.151 

2.182e-58 

33.224 

77.716 

814.438 

9.4 

4.329 

20.276 

2.023e-58 

33.272 

77.908 

818.332 

9.DL.1 

60.213 

97.206 

6.975e-56 

62.861 

196.261 

3214.621 

^774- 

9.1 

2.895 

8.239 

1.689e-83 

113.013 

103.383 

197.908 

9.2 

2.866 

9.236 

1.461e-83 

113.398 

104.921 

229.062 

9.3 

3.962 

10.424 

1.382e-83 

113.854 

106.748 

266.057 

9.4 

4.986 

10.552 

1.572e-83 

113.902 

106.941 

269.952 

9.DL.1 

13.807 

27.932 

5.415e-81 

106.562 

122.140 

819.291 

D 

wtm 

9.1 

1.873 

17.960 

2.287e-58 

32.757 

75.532 

759.839 

1 

9.2 

2.755 

18.959 

1.983e-58 

33.148 

77.099 

791.565 

■ 

9.3 

3.946 

20.145 

2.017e-58 

33.613 

78.959 

829.232 

l 

9.4 

4.328 

20.269 

1.871e-58 

33.662 

79.153 

833.166 

■ 

1 

9.DL.1 

60.213 

97.200 

6.533e-56 

63.789 

199.657 

3273.031 

I 

9.1 

1.613 

6.873 

1.416e-81 

50.935 

51.013 

195.564 

I 

9.2 

2.714 

7.971 

1.293e-81 

51.366 

52.735 

230.433 

9.3 

3.927 

9.155 

1.617e-81 

51.827 

54.582 

267.838 

1 

9.4 

3.636 

8.706 

1.345e-81 

51.653 

53.884 

253.698 

I 

9.DL.1 

14.948 

27.989 

6.016e-79 

59.203 

84.086 

865.279 

B 

9.1 

1.662 

6.327 

5.082e-84 

37.336 

40.060 

87.841 

B 

9.2 

2.717 

7.409 

4.386e-84 

37.761 

41.758 

222.226 

B 

9.3 

3.929 

8.588 

4.149e-84 

38.223 

43.608 

259.672 

9.4 

3.777 

8.234 

4.732e-84 

38.083 

43.047 

248.312 

9.DL.1 

25.998 

1.784e-81 

43.582 

69.630 

646.342 

2 

1 

9.1 

1.873 

17.957 

2.191e-58 

33.150 

76.730 

773.391 

9.2 

2.755 

18.954 

1.850e-58 

33.549 

78.324 

805.688 

9.3 

3.946 

20.141 

1.906e-58 

34.022 

80.218 

844.028 

9.4 

4.327 

20.264 

1.769e-58 

34.071 

80.414 

848.004 

9.DL.1 

60.213 

97.193 

6.096e-56 

64.737 

203.070 

3331.446 

7k 

9.1 

1.613 

6.318 

7.088e-82 

28.283 

32.813 

197.462 

9.2 

2.714 

7.416 

6.475e-82 

28.721 

34.566 

232.953 

9.3 

3.927 

8.599 

8.089e-82 

29.191 

36.446 

271.027 

9.4 

3.636 

8.151 

6.732e-82 

29.013 

35.736 

256.634 

9.DL.1 

14.948 

27.435 

3.015e-79 

36.699 

66.476 

879.136 

27/4 

9.1 

1.620 

5.990 

2.987e-84 

24.049 

29.071 

188.742 

9.2 

2.714 

7.086 

2.578e-84 

24.487 

30.823 

224.219 

9.3 

3.927 

8.264 

2.429e-84 

24.957 

32.703 

262.294 

9.4 

3.653 

7.826 

2.792e-84 

24.780 

31.996 

247.981 

9.DL.1 

13.805 

25.617 

1.069e-81 

31.334 

59.890 

822.696 

References 


Box,  G.  E.  P.  and  Wilson,  K.  B.  (1951).  On  the  experimental  attainment  of  optimum 
conditions.  Journal  of  the  Royal  Statistical  Society ,  Ser.  B,  13,  1—45. 

Box,  G.  E.  P.  and  Hunter,  J.  S.  (1957).  Multi-factor  experimental  designs  for  exploring 
response  surfaces.  Annals  of  Mathematical  Statistics,  28,  195—241. 

Box,  G.  E.  P.  and  Draper,  N.  R.  (1987).  Empirical  model-building  and  response  surfaces. 
John  Wiley  &  Sons,  New  York. 

Draper,  N.  R.  (1985).  Small  composite  designs.  Technometrics,  27,  173—180. 

Draper,  N.  R.  and  Lin,  D.  K.  J.  (1990).  Small  response-surface  designs.  Technometrics, 
32,  187-194. 

Ghosh,  S.  (1979).  On  robustness  of  designs  against  incomplete  data.  Sankhya  40,  Ser.  B, 
Pts  3  and  4,  204-208. 

Hartley,  H.  O.  (1959).  Smallest  composite  designs  for  quadratic  response  surfaces.  Bio¬ 
metrics,  15,  611—624. 

Khuri,  A.  I.  and  Cornell,  J.  A.  (1987).  Response  surfaces,  designs  and  analyses.  Marcel 
Dekker.  Inc.,  New  York. 

Plackett,  R.  L.  and  Burman,  J.  P.  (1946).  The  design  of  optimum  multifactorial  experiments. 
Biometrika,  33,  305—325. 

Westlake,  W.  J.  (1965).  Composite  designs  based  on  irregular  fractions  of  factorials.  Bio¬ 
metrics,  21,  324—336. 


RPR-04-1995  09=26 


DTIC 


703  274  9270  P.03/04 


PI  Name  (Last  First  MI)  Basu,  Asit  P. 

Institution  University  of  Missouri-Columbia 

Contract/Grant  No.  F49620-92-J-0371 

A.  Publications  in  peer-reviewed  professional  journals  and  refereed  book  chapters  during  the 
reporting  period,  1  October  1992  through  30  September  1993. 

Book  Edited: 

Basu,  Asit  P.,  Editor, 

ADVANCES  IN  RELIABILITY,  North-Holland,  Amsterdam,  1993. 

Book  Chapters: 

1 .  Klein,  John  P.  and  Goel,  Prem  K.,  Editors,  Survival  Analysis:  State  of  the  Art,  Kluwer 
Academic  Publishers,  Dordrecht,  1992. 

•Basu,  A.  P.,  Life  Testing  and  Reliability  Estimation  with  Asymmetric  Loss  Function. 

2.  Sen,  Pranab  K.  and  Salama,  Ibrahim  A.,  Editors,  Order  Statistics  and  Nonparametrics: 
Theory  and  Applications. 

•Basu,  A.  P.,  On  a  test  for  exponentiality  against  Monotone  Failure  Rates. 

3.  Goel,  Prem,  K.  and  Iyengar,  N.  S.,  Editors,  Bayesian  Analysis  in  Statistics  and 
Econometrics,  Springer- Verlag,  New  York,  1992. 

•Basu,  A.  P.,  Bayesian  Approach  to  Some  Problems  in  Life  Testing  and  Reliability 
Estimation. 

4.  Basu,  A.  P.,  Editor,  Advances  in  Reliability,  North-Holland,  1993. 

•Basu,  A.  P.,  Characterizations  of  a  Family  of  Bivariate  Exponential  Distributions. 

5.  Basu,  A.  P.,  Editor,  Advances  in  Reliability,  North-Holland,  1993. 

•Basu,  A.  P.,  Bayesian  Reliability  of  Stress-Strength  Systems. 

6.  Cothcm,  Richard  C.  and  Ross,  Phillip  N.,  Editors,  Environmental  Statistics  Assessment 
and  Forecasting. 

•  Basu,  A.  P.,  Some  Problems  of  “Safe  Dose”  Estimation. 


1 


APR-04-1995  09=26 

k  , 


DTIC 


703  274  9270 


B .  Number  of  additional  researchers  working  with  the  Principal  Investigator 

Faculty  one,  A.  P.  Basu 

Postdocs 
Graduate  Students 
Other 

C.  Professional  honors  received  by  all  of  the  above  listed  contributors 
Asit  P.  Basu  was  elected  Fellow  of  the  following  societies: 

1..  American  Association  for  the  Advancement  of  Science,  1987 

2.  American  Statistical  Association,  1983 

3.  Institute  of  Mathematical  Statistics,  1983 

4.  Royal  Statistical  Society,  England,  1974 
and  also  was  elected 

5.  Member  of  International  Statistical  Institute.  1987. 


2 


P . 04/04 


TOTAL  P.04 


